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The gauge-invariant generation of a dynamical, momentum-dependent gluon mass is intimately 
connected with the presence of non-perturbative massless poles in the vertices of the theory, which 
trigger the well-known Schwinger mechanism. In the deep infrared the integral equation that 
governs this effective gluon mass assumes a particularly simple form, which may be derived 
following two seemingly different, but ultimately equivalent procedures. In particular, it may be 
obtained either as a deviation from a special identity that enforces the masslessness of the gluon 
in the absence of massless poles, or as a direct consequence of the appearance of a non-vanishing 
bound-state wave function, associated with the details of the actual formation of these massless 
poles. In this presentation we demonstrate that, due to profound relations between the various 
ingredients, the two versions of the gluon mass equation are in fact absolutely identical. 
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1. Introduction 

It is well-established by now that the dynamical generation of an effective gluon mass [p ex- 
plains in a natural and self-consistent way the infrai^ed finiteness of the (Landau gauge) gluon 
propagator and ghost dressing function, observed in large volume lattice simulations for both 
SU (2) and SU (3) 0, |5|] gauge groups (for an alternative approach see [^]). Given the non- 
perturbative nature of the mass generation mechanism, the usual starting point in the continuum 
is the Schwinger-Dyson equations (SDEs) governing the Green's functions under scrutiny. In the 
framework provided by the synthesis of the pinch technique (PT) [JT), ^, |8]| with the background field 
method (BFM) [Q], these complicated integral equations are endowed with a variety of important 
properties, which allow a much tighter control on the truncations adopted and the approximation 
schemes employed. 

Probably the most crucial theoretical ingredient in this context is a special type of vertex, 
denoted by V , which contains massless, longitudinally coupled poles. This vertex complements the 
all-order three-gluon vertex entering into the SDE governing the gluon self-energy, and captures 
the underlying mass generating mechanism, which is none other than the Schwinger mechanism. 
Specifically, the QCD dynamics is assumed to generate massless bound-state excitations, which, in 
turn, give rise to the aforementioned poles that appear inside the vertex V . 

At the level of the SDE for the gluon propagator, the analysis finally boils down to the deriva- 
tion of an integral equation that governs the evolution of the dynamical gluon mass, m^{q^), as a 
function of the momentum in a way similar to the more familiar case of the dynamical gen- 
eration of a constituent quark mass. However, unlike what happens with the quark gap equation, 
which, due to its Dirac structure, is unambiguously separated into two equations, governing the 
quark mass and wave-function, the derivation of the corresponding equations for the gluon mass 



and wave-function requires the introduction of some additional kinematic criteria [|10|]. For the pur- 
poses of this presentation, we will focus on the integral equation for m^{q^) in the infrared limit, 
i.e, as — )• 0, where this separation becomes unique and totally unambiguous: one assigns to the 
mass equation all contributions that do not vanish as q^ 0. 

It turns out that the equation for m^(0) may be derived following two rather distinct procedures, 
which eventually express the answer in terms of seemingly different quantities. Roughly speaking, 
the first procedure, which is operationally closer to the standard SDE treatment, identifies the con- 
tribution to the gluon mass equation as the deviation produced to the so-called "seagull-identity" by 
the inclusion of the vertex V . Interestingly enough, this can be accomplished without knowledge of 
the explicit closed form of the vertex V , relying only on its general features, most notably its lon- 
gitudinal nature and the Slavnov-Taylor identities that it satisfies. The second procedure expresses 
the answer in terms of quantities such as "the bound-state wave-function", denoted by B\, appear- 
ing in the study of the Schwinger mechanism from the point of view of the actual formation of the 
required bound state excitations. This particular quantity satisfies a homogeneous Bethe-Salpeter 
equation, which, as has been shown recently, admits indeed non-trivial solutions [[III]. 

To be sure, these two procedures, despite their apparent differences, must be ultimately equiv- 
alent. In this talk we will argue that this is indeed the case, by showing that the two mass equations 
obtained coincides in the deep infrared limit, thus providing an important self-consistency check 
for this entire approach. 
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Figure 1: The "one-loop dressed" gluon contribution to the PT-BFM gluon self-energy. White (black) 
circles denote fully dressed propagators (vertices); a gray circle attached to the external legs indicates that 
they are background gluons. Within the PT-BFM framework these two diagrams constitute a transverse 
subset of the full gluon SDE. 



2. Dynamical gluon mass: general concepts and ingredients 

Let us start by setting up the notation and reviewing some of tlie most salient features of tlie dy- 
namical gluon mass generation scenario formulated witliin the PT-BFM framework. The full gluon 
propagator A'^^^,(^) = 5"^A^v('?) in the Landau gauge is defined as 

A^v(^) = -iPi,v{q)Kq^); P^viq) = g^v - (2.1) 

where the scalar factor A(<7^) and the gluon self-energy n^v(^) are related through 

A-\q^)=q^ + iniq^y, U^yiq) = P^y{q)n{q^). (2.2) 

For later convenience let us also define the inverse gluon dressing function, J{q^), as 

K-\q^)=q^J{q^); (2.3) 

then, in the presence of a dynamically generated mass, a massive gluon propagator is naturally 
described in term the natural form of the expression 

KXq^) = q^ln{q^)-m\q^): (2-4) 

where the first term corresponds to the "kinetic term", or "wave function" contribution, whereas 
the second is the (positive-definite) momentum-dependent mass. Notice that the symbol 7,„ in- 
dicates that effectively one has now a mass inside the corresponding expressions: for example, 
whereas perturbatively J{q^) ~ In^^, after dynamical gluon mass generation has taken place, one 
has J,n{q^) ~ ln(^^ + m^). 

The usual starting point of our dynamical analysis is the SDE governing the gluon propagator. 
Within the PT-BFM framework that we employ, one may safely truncate the SDE series down to its 
"one-loop dressed version" containing gluonic contributions only, given by the diagrams (ai ) and 
(02) shown in Fig. [l] [|l2[, |l3[ [T^. In fact, the resulting (approximate) gluon self-energy n^v(^) is 



transverse, i.e., 

qf'n^y{q)=0. (2.5) 
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The PT-BFM equation for the conventional propagator reads, in this case, 



A-\q')P^'{q) 
where, in the Landau gauge. 



[l+G(^2)]2 



(2.6) 



A',\q) 



A^i^) 



P"^{k) 



(k + qy 



A{k), 



P''^k)-^-^--^A{k), 



{k + qf 



d 

k^{k + qy 
k^ik + qf 



A{k), 



(2.7) 



and the J-dimensional integral measure (in dimensional regularization) is defined as = J d^k. 
The vertex IT is the fully-dressed BQQ vertex, connecting a background gluon (B) to two quantum 
gluons (Q), which naturally appears in the PT-BFM framework and that has been studied in detail 
in [ jlsl ] (see also below); in addition. 



^fap P) = Sap {r-p)^+ gpti {p-q)a+ gafi ( 



(2.8) 



Finally, the function G{q^) represents the scalar co-factor of the g^y component of the special 
two-point function A^y{q), defined as 



^nv{q) = -ig^CA A^{k)D{q-k)Hya{-q,q-k,k) 

Jk 



g,Mq') + '-^Liq 

q^ 



(2.9) 



where we have introduced the ghost propagator D"^{q^) = 8'^^D{q^), which is related to the ghost 
dressing function F{q^) through D{q^) = ^-^^ Notice that in the Landau gauge, an important 
exact (all-order) relation exists, linking G{q^) and L{q^) to the ghost dressing function F{q^), 
namely [16, 17] 

(2.10) 



F-\q') = \+G{q')+L{q^); 



in addition G coincides with the well-known Kugo-Ojima function []16[ |17p . 

The Schwinger mechanism [18] is integrated into the SDE of the gluon propagator through 
the form of the three-gluon vertex [19, 20, |^]. In particular, a crucial condition for the real- 
ization of the gluon mass generation scenario is the existence of a special vertex, to be denoted by 
Va^v{q,r,p), which must be completely longitudinally coupled, i.e., must satisfy 



P"'-{q)P^'^ir)P^'^ip)Va,Aq,r,p)=0. 



(2.11) 
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Figure 2: The W three-gluon vertex. Thick (red onUne) internal gluon lines indicates gluon propagators 
with an effective mass. 

This vertex is instrumental for maintaining gauge invariance, given that the massless poles that it 
must contain in order to trigger the Schwinger mechanism, act, at the same time, as composite, 
longitudinally coupled Nambu-Goldstone bosons. Specifically, in order to preserve the gauge in- 
variance of the theory in the presence of masses, the vertex Va^vi^i ^iP) must be added to the BQQ 
(fully-dressed) three-gluon vertex ^a^ivi^l^^'^p)^ giving rise to the new full vertex, Wfy,^^{q,r,p), 
defined as [0] 

^'a^viq^r.p) = Jra^v{q,r,p)+Va^v{q,r,p). (2.12) 

To see in detail how gauge invariance is preserved, notice that when the Schwinger mechanism is 
turned off, the BQQ vertex alone satisfies the WI 

q"Ta^yiq,r,p)=p^Jip^)P^y{p)-r^Jir^)P^y{r), (2.13) 

when contracted with respect to the momentum of the background gluon. By requiring that 

q"Va^,viq,r,p)=mHr^)P^y{r)-m\p^)P^^{p), (2.14) 

we see that, after turning the Schwinger mechanism on, the corresponding WI satisfied by IT' would 
read 

(l"^anviq,r,p) = q"[jr{q,r,p)+V{q,r,p)]^^^ 

= [p^J{p^)-m\p^)]P,,{p) - -m\?)]P,,{r) 

= A^\p^)P,,{p)-A;„\r^)P^,{r), (2.15) 

which is indeed the identity ( ^.13| ) satisfied by IT with the replacement of the conventional gluon 
propagators appearing on the rhs by massive propagators: A^' — )• The remaining (more diffi- 
cult) STIs, triggered when contracting Jr'g,^y{q,r,p) with respect to the other two legs are realized 
in exactly the same fashion ||T^. 

The next step is to insert IT'^^^, r, p) into the SDE equation satisfied by the gluon propagator, 
thus obtaining the new SDE of Fig. ^ From the resulting SDE one can obtain two separate equa- 
tions, the first one governing the behavior of J,„{q^), and the second one describing the dynamical 
mass np-{q^). The general idea is the following: the terms appearing on the rhs of the SDE may be 
separated systematically into two contributions, one that vanishes as <7 — )• and one that does not; 



5 



Infrared properties of the gluon mass equation 



Daniele Binosi 







Figure 3: Diagrammatic representation of the gluon one-loop dressed diagrams before and after dynamical 
gluon mass generation has taken place: the propagators and vertices on the rhs have now become A„, and IT'. 



the latter contribution must be set equal to the corresponding non- vanishing term on the Ihs, namely 
—m^{q), while the former will be set equal to the vanishing term of the Ihs, namely q^Jmiq^)- 

The exact way how to implement this particular separation is rather elaborate, and requires the 
use of the so-called "seagull identity" [ ^ ] 

d 



[k'A'{k) + ^ [A{k)=0, 
Jk 2 Jk 



(2.16) 



where the "prime" denotes differentiation with respect to k^. This identity plays a a crucial role, 
forcing all possible quadratic divergences, associated with integrals of the type /^A(^), and varia- 
tions thereof, to cancel exactly. In the case of the dimensional regularization that is used through- 
out, the presence of such integrals would give rise to divergences of the type mQ^l/e), where mo 
is the value of the dynamically generated gluon mass at = 0, i.e., niQ = m(0); if a hard cutoff A 
were to be employed, these latter terms would diverge quadratically, as A^. The disposal of such 
divergences would require the introduction in the original Lagrangian of a counter-term of the form 
wIqA^, which is, however, forbidden by the local gauge invariance, which must remain intact. No- 
tice that the two types of integrals appearing on the Ihs of Eq. ( 2.16| ) are individually non- vanishing 
(in fact, they both diverge); it is only when they come in the particular combination shown above 
that they sum up to zero. 



The seagull cancellation implemented by Eq. ( 2.16 ) is a point of paramount importance, be- 
cause it further re-enforces the fundamental assertion that the gluon mass generation, when im- 
plemented correctly, is absolutely compatible with the underlying BRST symmetry. Indeed, in the 
mass generation picture advocated in a series of recent articles [ 14 , 25 ] the Lagrangian of 

the Yang-Mills theory (or that of QCD) is never altered; therefore, the only other possible ways of 
violating the gauge (or BRST) symmetry would be (/) by not respecting, at some intermediate step, 
some of the Wis and STIs satisfied by the Green's functions involved; for example, in the conven- 
tional SDE formulation, a naive truncation would compromise the transversality of the resulting 



gluon self-energy, i.e., the fundamental transversality condition ( |2.5[ ) would be no longer valid; 
and/or (//) by introducing seagull-type divergences, of the type mentioned above. Evidently, nei- 
ther (/) nor (//) happens, thanks to the powerful field-theoretic properties encoded into the PT-BFM 
formalism that we employ. 

Returning to the gluon SDE, a rather elaborate analysis [ p^ ] gives rise to two coupled integral 
equations, one for Jm{q) and one for the momentum-dependent gluon mass, of the generic type 



m 



(q^) = 1+ J^Jfi{q\m\A,„), 
(q^) = J^je2{q^,m\A„,). 



(2.17) 
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such that q^Ji^i {q^ , , A^) — >• 0, as — > 0, whereas J^2 {q^ , , A„, ) 7^ in the same lunit, precisely 
because it includes the term \/q^ contained inside Va^v{Q,r,p). 

There is a relatively straightforward way to derive the closed form of the gluon mass equation 
in the limit q^ — )• 0, without knowledge of the explicit form of V; all one needs is to postulate its 
existence and the properties mentioned above, most notably the identities ( 2.11 ) and ( |2.14 ). In 
order to do that, we start by noticing that, since the vertex IT' satisfies the identity ( 2.15 ), the gluon 
self-energy is transverse even in the presence of masses; as a result, restoring the Lorentz structure 
into the second equation of ( [2.17[ ), we have that 

m\q^)P^y{q) = P^y{q) J^je2{q\m\A^). (2.18) 

Now, due to its longitudinal nature, the vertex V can only furnish the part on the rhs proportional 
to q^qy/q^. The question is, where will the g^y part come from. A bit of thought reveals that this 
term can only emerge as a deviation from the seagull identity, which enforces the masslessness of 
the gluon when V = 0. In fact, this identity, given its nature and function, can only operate among 
terms that are proportional to g^y; this is so because the basic seagull contribution proportional to 
/j(.A(^) stems entirely from graph {a2), which has no momentum dependence, i.e., it can only be 
proportional io g^y. 

Specifically, the required contribution proportional to stems from the term 

/ ; Mi\Mi ^ Ak + qfJ,n{k + q)-k^Jm{k) ^ . 2n , ^M^v „ . 2n .oiox 

k^kyA{k)A{k + q) =8fivCi{q ) + -^C2{q ), (2.19) 

which originates from A^^(^). Then, since for any function f{k^) one has the result 

[cos^df{k^) = ^ [f{k^), (2.20) 
Jk a Jk 

it is easy to demonstrate that 

Ci(0) = J^k^A^{k)[k^J{k)]' 

= jk^A^{k)[A-\k)+m^{k)\ 

= f k^A^{k)[m^{k)]' - [ k^A'{k). (2.21) 
Jk Jk 

The first term on the rhs forms part of g^yJ^{q-^,m^,A„,), while the second, after taking all relevant 
multiplicative factors correctly into account, cancels against the contribution of graph (ai), by 
virtue of the seagull identity (2.16). This basic observations can be generalized to include the 
additional structures coming from the remaining terms in (2.7), and in particular A'^^{q), whose 
net effect is the one of reducing by one power the 1 + G factor appearing in the denominator of 
Eq. (|!|) lH^. 

Then, restoring all relevant factors, and using the fact that in 4 dimension L(0) = [ [l7| ] so that 
by virtue of Eq. ( 2.10 ) one can trade the 1 + G combination for the inverse of the ghost dressing 
function, one finally obtains, in Euclidean space, 

m^O) = lg^CAF{0) [k^[m\k)]'A^{k) 
2 Jk 

= -'ig^CAF{0) jn?{k)A{k)[k^A{k)\' . (2.22) 
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Figure 4: The kernel [q^A^{q^)]' obtained from the SU{3) (left) and SU{2) (right) lattice data. In the 
inset we show the corresponding lattice results for the gluon propagator, renormahzed at /i = 4.3 GeV and 
jU — 2.2 GeV respectively. 



and, after carrying out the angular integration, the final equation 

m\0) = -^as.f (0) ^Aym\y)[y^^\y)t (2.23) 

where = /An and y = k^ (and therefore the "prime" indicates now derivatives with respect 
toy). 

It turns out that the specific form of the mass equation ( 2.23 ) introduces a non-trivial constraint 
on the precise behavior that A must display in the region between (1-5) GeV^. Specifically, in order 
for the gluon mass to be positive definite, the first derivative of the quantity q^A{q^) (the "gluon 
dressing function") must furnish a sufficiently negative contribution in the aforementioned range 
of momenta. Interestingly enough, the A obtained from the lattice, shown in Fig. ^ has indeed this 
particular property; in the plots shown in Fig. ^ we clearly see that both derivatives change their sign 
in the intermediate momenta region, which constitutes precisely the required behavior. This is to be 
contrasted to what happens, for example, in the case of a simple massive propagator l/{q^ + m^) 
or with the Gribov-Zwanziger propagator jicf" + m'*) (with m constant); the derivatives of the 
corresponding dressing functions, j {q^ + m^) and cf jicf' + m'^), respectively, are positive in the 
entire range of (Euclidean) momenta, thus excluding, in this context, the possibility of a positive- 
definite gluon mass. 



3. Schwinger mechanism in the bound-state language: a self-consistency check 



The key assumption when invoking the Schwinger mechanism in Yang-Mills theories, such 
as QCD, is that the required poles may be produced due to purely dynamical reasons; specifically, 
one assumes that, for sufficiently strong binding, the mass of the appropriate bound state may be 



reduced to zero [20| , [21| , [22| ]. It is precisely these poles that constitute the main ingredient of 
the special vertex V. Specifically, all terms appearing in this vertex are proportional to l/q^, l/r^, 
l/p^, and products thereof, and can be separated in two distinct parts, according to [ITl]] 



Vanv{q,r,p) = Uayiv{q,r,p) +Ranv{q,r,p) 



(3.1) 



8 



Infrared properties of the gluon mass equation 



Daniele Binosi 




m, ^ 



Figure 5: (A) The vertex UapiV is composed of three main ingredients: the transition ampHtude, la, which 
mixes the gluon with a massless excitation, the propagator of the massless excitation, and the massless 
excitation-gluon-gluon vertex. {B) The Feynman rules (with color factors included) for (/) the propagator 
of the massless excitation and (n) the "proper vertex function", or, "bound-state wave function", Bjiv. 



The R term contains structures proportional to r^/r^ and/or Pv/p^, which vanish when contracted 
with the product Pn'^{r)Pyiy{p) coming from the Landau gauge gluon propagators, and therefore 
are not relevant for us. On the other hand, the U term can be written as [see Fig. ^ (A)] 

Uativi(],r,p) = ila{q) B^y{q,r,p), (3.2) 

where the factor // represents the propagator of the scalar massless excitation, while the nonper- 
turbative quantity 

B^v{q,r,p) = Big^v + B2q^qv + Bjp^Pv + B4r^qv + Bsr^pv, (3.3) 

is the effective vertex describing the interaction between the massless excitation and two gluons 
[Fig. ^ (B)]. Finally, lai^l) = qaK^) the (nonperturbative) transition amplitude introduced in 
Fig. [s], allowing the mixing between a gluon and the massless excitation' . li is important to notice 
that, due to Bose symmetry with respect to the quantum legs, i.e., pL -i^ v and p r one has 
B\;i{q-,r,p) = — Bi 2(<7,p,?'), and therefore these two form factors vanish in the limit, g — 0. 

The WI ( 2.14| ) furnishes an exact relation between the dynamical gluon mass, the transition 



amplitude at zero momentum transfer, and the form factor B\. Specifically, contracting both sides 
of the WI with two transverse projectors, one obtains, 

P^'^{r)P'''''{p)q^Va^,{q,r,p) = [m\r) - m\p)]P^ {r)P^''' {p) . (3.4) 



' Notice that we do not absorb the extra factor of "i" coming from the Feynman rule in Fig.^into the definition of the 
transition function , as was done in | [n| ] ; this definition has the advantage that the function / is real in the Euclidean 
space. 
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Figure 6: The "squared" diagram. 



On the other hand, contracting the full expansion of the vertex ( ^ ) by the same transverse projec- 
tors and next contracting the resulting expression with the momentum of the background leg, we 
get 

q"P^''^{r)P''\p)Va^,{q,r,p) = -l{q)[B,g^, +B2q^q,]P^'^r)P''\p), (3.5) 
Thus, equating both results, one arrives at the exact relations 



I{q)Bi{q,r,p) = {p) - n? {r); B2{q,r,p) = 0. 



(3.6) 



Finally, carrying out the Taylor expansion around ^ = of both sides of the first relation in (p.q), 
and using the fact that the form factors B\ 2 vanish in this limit, we arrive at (Minkowski space) 



[m\p)r = -mB[{p), 



where we have defined 



B[ip)=B[{0,-p,p) = lim 



dBi{q,-p-q,p) 

dip+qy 



(3.7) 



(3.8) 



and have implicitly assumed that 7(0) is finite. 

We will next approximate the transition amplitude Ia{q), connecting the gluon with the mass- 
less excitation, by considering only diagram {d\) in Fig. ^, corresponding to the gluonic "one-loop 
dressed" approximation; we will denote the resulting expression by Ia{l)- In the Landau gauge, 
the amplitude /« (q) reads 

Ia{q) = hcAfmA{k + q)rfp^P^^{k)pP''{k + q)B^y{-q,-k,k + q^ (3.9) 



where the origin of the 1 /2 factor is combinatoric, and r|^^^ is the tree-level vertex of Eq. (2.8). 

At this point we are in the position to determine how the transition function is related to 
the dynamically generated gluon mass which originates from the inclusion of the vertex V in the 
corresponding gluon SDE. Indeed, at the level of approximation employed here we find that the 
contribution to the gluon self energy is 



Ar(^)lv = \g'CAfAiq + k)A{k)r^^fp-P{k)P^<'{k + q)U;, 
\ca jA{q + k)A{k)Tf^^P^P {k)pP- {k + q)B 



iv{-q) 



■q^l 272/ ^ 

r 



(3.10) 
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where we have used Eq. (3.2) (with ly ly), together with the property ly{—q) = —Iv{q) as well 



as Eq. ( |3.9| ). Notice that as anticipated this contribution is purely longitudinal. 

Since the inclusion of V in A3 has the same effect previously found, namely to effectively 
remove one power of ( 1 + G) in the denominator^ , we obtain the positive-definite result 



m\0)=g^FiO)PiO). 



(3.11) 



Of course, for the PT-BFM framework to be self-consistent, the two infrared equations for the 
gluon mass ought to coincide. To see that this is indeed the case, consider again the original mass 



equation (2.22) , and let's insert the relation ip7T\) - at the same level of approximation used in the 



derivation of (3.11), i.e., with the replacement 7(0) — > 7(0) - into the integral in the first line, to get 

]g^CAF{0)Je[mHk)]'A\k) 



gV(0)7(0) 



On the other hand, a a bit of algebra reveals that, in Euclidean space, on has [ ]llp 

3 



7"(0) = -/Ca [ k^A^{k)B\{k), 
2 Jk 



so that we indeed ( |3.12[ ) gives rise to the relation 

mHO)=g^FiO)PiO), 



(3.12) 



(3.13) 



(3.14) 



which is none other than Eq. (3.11). 



4. Conclusions 

In this work we have presented a highly non-trivial self-consistency check related to the integral 
equation describing the dynamics of the gluon mass in the deep infrared. 



We started by showing how relation (2.22) can be derived only by assuming the existence of 



the special vertex V, without providing any information either on its own diagrammatic structure, 
nor on its composition in terms of other dynamical quantities, such as the bound-state wave function 
B. Next, we presented an alternative formulation, where all the intricate dynamical ingredients that 
trigger the Schwinger mechanism, subject to the constraints imposed by gauge invariance, were 
furnished. Both methods give rise to the same equation, proving that they represent nothing but 
different facets of the same underlying dynamics. 

In the future it would be particularly important to generalize the consistency check described 
to the full momentum range, i.e., extend it to the full dynamical equation governing the evolution 
of the mass nP'{q^) and not only to the deep infrared limit m^(0). In addition, one should go beyond 
the one-loop approximation by including the missing two-loop dressed gluon diagrams in both the 
SDE as well as the bound-state analysis. Work along these directions is already in progress. 



^It is precisely the inclusion of this term that account for the difference of Eq. ( 3.1 1 ) with respect to the analogous 
equation (3.38) of^. 
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